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ДМИТРИЙ РЕЧКИН 
СОЛНЕЧНАЯ СИСТЕМА И СОЛНЕЧНАЯ АКТИВНОСТЬ 


1. МАТЕРИАЛ И МЕТОДЫ 


Гипотезы о воздействии планет Солнечной системы на проявления солнечной активности вы- 
двигались неоднократно!?’3”4, при этом вопрос о физической природе процессов, определяю- 
щих наблюдаемые проявления солнечной активности - либо не ставился вообще, либо решал- 
ся с использованием частных, ограниченных моделей и подходов. 


В рамках данной работы предлагается следующий подхол: 


Солнце, как звезда спектрального класса С2\/, характеризуется наличием внутреннего 
ядра, окружающей его зоны лучистого переноса, переходящей в конвективную зону и 
атмосферу. В составе последней выделяются фотосфера, хромосфера, корона и сол- 
нечный ветер. Доля водорода в массе Солнца составляете 73,46%. Несмотря на доста- 
точно высокую (1,57.10’К) температуру солнечного ядра, возможно, что при сверхвысо- 
ких давлениях его электрическая проводимость близка к высокотемпературной сверх- 
проводимости. 


Сверхпроводимость и сверхтекучесть взаимосвязаны”, поверхность (граница раздела 
фаз) солнечного ядра должна демонстрировать неустойчивость в гидродинамическом 
смысле. Поэтому резонно ожидать, что эта поверхность весьма подвержена деформа- 
циям в результате различных флуктуаций (в частности, гравитационных). 


Оценить периодичность и масштаб флуктуаций возможно, поскольку они ведут к «отры- 
ву» от ядра высокоэнергетических «сгустков», постепенно поднимающихся от ядра к фо- 
тосфере; по мере продвижения от ядра к поверхности звезды эти сгустки принимают то- 
роидальную форму, как наиболее оптимальную при движении в текучей среде. 


На протяжении всего пути до поверхности тороидальная форма сгустков сохраняется, а 
за счет высоких внутренних токов поддерживается их более высокая (по сравнению с 
окружающим веществом) энергонасышенность. Достигнув фотосферы, тороидальные 
сгустки прорывают её, образуя видимые солнечные пятна (существующие, как уже дока- 
зано, в виде парных структур). 
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Таким образом, чем больше флуктуаций, деформирующих ядро Солнца -— тем больше всплы- 
вающих сгустков и, соответственно, больше пятен. Иными словами, исследуя динамику появ- 
ления пятен, можно составить представление о динамике сил, влияющих на энергетические 
выбросы со стороны солнечного ядра (и, как несложно догадаться, опосредованно - на общие 
показатели динамики солнечной активности). 


В работе использованы методы многомерной статистики, в частности — метод главных компо- 
нент8’’10, Все необходимые расчёты выполнены с применением калькулятора многомерной 
статистики $ТАТ!С/ЕХ!11 собственной разработки. Программы обработки данных на языке 
ЯКОБИ!? приводятся в приложении, исходные данные доступны по указанным ссылкам. 


1.1. ПокаАЗАТЕЛИ СОЛНЕЧНОЙ АКТИВНОСТИ 


Использовались исходные данные среднегодовых и среднемесячных значений чисел Вольфа, 
опубликованные на сайте Мирового центра данных $15013. В соответствии с представлениями 
о смене магнитной полярности циклов солнечной активности было дополнительно проведено 
простое преобразование для представления показателей активности с учётом полярности цик- 
ла (поочередная смена знака чисел Вольфа, по циклам). 


Дополнительно вычислялись (как описано в предыдущей работе“) «Фактор Е» (наблюдаемый 
размах активности) и «Фактор Х» (скрытая энергетика активности), как с учётом полярности 
цикла, так и без её учёта. 


1.2. ПАРАМЕТРЫ ДВИЖЕНИЯ СОЛНЦА 


Основываясь на взглядах Теодора Ландшейлдта!?'16, по данным эфемерид Солнца и планет Сол- 
нечной системы, Карл Смит (Са ЗтКВ) рассчитал и опубликовал таблицы!7 данных по смеще- 
ниям барицентра Солнечной системы относительно центра Солнца. Таблицы включают данные 
(за период с 3000 г. до н.э. по 3000 г. н.э.; всего за 6000 лет, с интервалом 5 суток между отсчё- 
тами): 


Т, дата отсчёта (год в десятичном представлении); 

А; склонение центра Солнца относительно барицентра Солнечной системы (в градусах); 

А; — отклонение центра Солнца от барицентра Солнечной системы, в долях радиуса Солнца 
(Вот = 6,957000-10:° см): 
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Г, момент импульса, определяется как [ = А-р = А-Млм, где А — отклонение центра Солнца от 
барицентра Солнечной системы, М - масса Солнца, у — скорость орбитального движения 
(Мл = 1,98847-103° г); 

т; крутящий момент, определяется какт = Р.А, где Е - действующая сила, К — отклонение 
центра Солнца от барицентра Солнечной системы, вычисляется как отношение прира- 
щения момента импульса за период наблюдений к длительности этого периода. 


Для расчётов в рамках настоящей работы использовались усредненные (сглаженные) показа- 
тели отклонения центра Солнца от барицентра (№), момента импульса (Ё) и крутящего момента 
(т), приведенные к среднемесячным и среднегодовым значениям (по исходным данным Карла 
Смита). 


2. РЕЗУЛЬТАТЫ 


2.1. АВТОКОРРЕЛЯЦИИ РЯДОВ ПОКАЗАТЕЛЕЙ ДВИЖЕНИЯ СОЛНЦА 


Автокорреляции (корреляции значений временного ряда с самим собой в предшествующие 
годы, в диапазоне лага лет от 0 до 300) на интервале 6000 лет представлены на графике (Рис.1). 
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Рис.1. Автокорреляции временных рядов - годичное сглаживание: 
В (черная линия и маркеры), Ё (синие), т (красные) 


Ко всем показателям с их временными отсчетами применялось годичное сглаживание. Для ря- 
да значений отклонения центра Солнца от барицентра Солнечной системы максимальное зна- 
чение автокорреляции составляет 0,991 при лаге до 180 лет; для ряда значений момента им- 
пульса максимальная автокорреляция составляет 0,994 при лаге до 180 лет; для ряда значений 
крутящего момента максимальная автокорреляция составляет 0,991 при том же лаге до 180 
лет!8. 


Для последующих расчетов важно определить средний период цикличного изменения показа- 
телей. Используя данные с годичным сглаживанием при лаге 180 лет, получены такие оценки: 


е длЯ ряда значений отклонения центра Солнца от барицентра Солнечной системы сред- 
ний период составил Тр = 19,73 года; 
® длЯ ряда значений момента импульса средний период Т, = 19,8 года; 


18 Программа сбип-уеаг.фас (см. Приложение П1) 


е для ряда значений крутящего момента Т, = 18,5 лет (по всей видимости, сказывается 
влияние нелинейного характера изменений этого параметра во времени). 


Аналогичные расчёты? проведены для показателей с месячным сглаживанием, однако в каче- 
стве максимального значения лага было выбрано 600 (месяцев, или 50 лет). Для ряда значений 
отклонения центра Солнца от барицентра Солнечной системы первое максимальное положи- 
тельное значение автокорреляции составило 0,426 при лаге 261 месяц (21,7 года). Для ряда 
значений момента импульса первое максимальное положительное значение автокорреляции 
составило 0,590 при лаге 244 месяца (20,3 года). Для ряда значений крутящего момента первое 
максимальное положительное значение автокорреляции составило 0,413 при лаге 273 месяца 
(22,7 лет). 
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Рис.2. Автокорреляции временных рядов - месячное сглаживание: 
В (черная линия), Ё (синяя), т (красная) 


Представленные графики отражают автокорреляции названных временных рядов с помесяч- 
ным сглаживанием, в диапазоне лага месяцев от 0 до 600 (от 0 до 50 лет). 


2.2. ЦИКЛИЧНОСТЬ РЯДОВ ПОКАЗАТЕЛЕЙ ДВИЖЕНИЯ СОЛНЦА 


В расчётах и демонстрации цикличности процессов, определяющих движение Солнца, исполь- 
зованы данные с годичным сглаживанием при лаге 180 лет. Обработка выполнена методом 
главных компонент по стандартной для временных рядов схеме. Построены фазовые портреты 
(проекции траектории процесса на фазовые плоскости главных компонент, в качестве объектов 
рассматривались годы). 
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Рис.3. Фазовые портреты ряда смещений центра Солнца относительно барицентра Солнечной системы: 
а - 1[-П, Б- Ш-ТУ, с -ТУ-И, а- УП-УШ главные компоненты 
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На Рис.3 показаны проекции траектории процесса годичных смещений центра Солнца относи- 
тельно барицентра Солнечной системы (за 3000 лет). С небольшими вариациями все цикличе- 
ские колебания представляют собой гармонические функции вида $т (С05), поскольку для них 
выполняется соотношение х2:+у2; = соп$ё, где х; и у; — координаты проекции точки (года) на оси 
соответствующих главных компонент. Для наглядности в каждой из проекций выделен цветом 
один цикл, для которого отмечены точки, соответствующие годам. Приходящиеся на |-\МШ глав- 
ные компоненты доли дисперсии, а также цикличность соответствующих временных рядов по- 
казателей движения Солнца приведены в Табл.1. 


На Рис.4 и Рис.5 показаны, соответственно, проекции траектории процессов годичных измене- 
ний момента импульса и крутящего момента движения Солнца. 
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Рис.4. Фазовые портреты ряда момента импульса движения Солнца: 
а - 1[-П, Б- Ш-ГУ, с -У-И, а- УП-УШ главные компоненты 
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Рис.5. Фазовые портреты ряда крутящего момента движения Солнца: 
а - 1-1, БВ - Ш-ГУ, с - У-ИГ а- УП-УШ главные компоненты 


Таблица 1 
Доли дисперсии главных компонент и цикличность соответствующих временных рядов 


Главные компоненты 


ии му 
Отклонение центра Солнца относитель- 33,50% | 33,04% | 9,96% | 9,83% | 2,78% | 2,76% 1,16% 


но барицентра Солнечной системы 
66,54% 19,79% 


19,8 года 
Момент импульса движения Солнца 


14,67% 


19,8 года 12,6 года 13,5 года 35,4 года 


Крутящий момент движения Солнца 31,21% | 30,77% | 14,03% | 13,89% | 3,95% | 3,93% | 0,35% 
61,98% 27,92% 7,88% 
19,8 года 12,6 года 13,5 года | 73а = 





Таким образом, при исследовании цикличности временных рядов параметров движения Солн- 
ца в Солнечной системе методом главных компонент выявлено, как минимум, четыре периода 
гармонической цикличности (по убыванию доли дисперсии): 19,8 года, 12,6 года, 13,5 года и 
35,4 года. Некоторый интерес может также представлять выявленная цикличность с периодом 
7,3 года, проявляющаяся при исследовании временного ряда значений крутящего момента 
движения Солнца. Ни один из этих периодов не совпадает с известными периодами обраще- 
ния отдельных планет вокруг центрального светила. Следует напомнить, что точность опреде- 
ления периодов невысока, поскольку в основе расчётов лежат временные ряды с годичным 
сглаживанием значений. 


2.3. ЗАПАЗДЫВАНИЕ В ПРОЯВЛЕНИЯХ СОЛНЕЧНОЙ АКТИВНОСТИ 


Перед выполнением финальных расчётов влияния движения Солнца на цикличность солнечной 
активности важно определить запаздывания (фазовые сдвиги) для отдельных временных рядов 
показателей, вычисленных на основе эфемерид (отклонения центра Солнца от барицентра 
Солнечной системы, момент импульса и крутящий момент). Поиск оптимальных значений за- 
паздывания (в годах) выполняется по максимальным значениям корреляции рядов показате- 
лей солнечной активности (числа Вольфа, как с учётом полярности цикла, так и без нее, факто- 
ры «Е» и «Х») с рядами параметров движения Солнца. Поскольку вычисления? проводятся для 
рядов в ограниченном интервале лет (1700 - 2019), максимальное значение лага приходится 
ограничить (до 50 лет). 


На Рис.6 показано изменение значений коэффициента корреляции между проявлениями сол- 
нечной активности и отклонением центра Солнца от барицентра Солнечной системы, в зависи- 
мости от лага (интервала запаздывания). 
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Рис.6. Корреляции проявлений солнечной активности 
с отклонениями центра Солнца от барицентра Солнечной системы 


Цвет линий обозначает: 
е зеленый — среднегодовые чисел Вольфа, без учёта полярности цикла; 
® красный — среднегодовые значения чисел Вольфа, с учётом полярности цикла; 


® синий — вычисленные значения «Фактора Е», без учёта полярности цикла; 
е чёрный — вычисленные значения «Фактора Х», без учёта полярности цикла. 
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На Рис.7 показано изменение значений коэффициента корреляции между проявлениями сол- 
нечной активности и моментом импульса движения Солнца, В зависимости от лага (интервала 
запаздывания), цветовые обозначения -— те же. 
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Рис.7. Корреляции проявлений солнечной активности с моментом импульса движения Солнца 


На Рис.8 показано изменение значений коэффициента корреляции между проявлениями сол- 
нечной активности и крутящим моментом движения Солнца, в зависимости от лага (интервала 
запаздывания), цветовые обозначения - те же. 
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Рис.8. Корреляции проявлений солнечной активности с крутящим моментом движения Солнца 


В качестве оптимальных значений лага (количества лет запаздывания) для дальнейших расчё- 
тов выбраны: 


для отклонения центра Солнца от барицентра Солнечной системы - 7 лет (по макси- 
мальным корреляциям со среднегодовыми числами Вольфа, с учётом полярности цикла; 
В7= 0,319) 


для момента импульса движения Солнца - 7 и 8 лет (по максимальным корреляциям со 
среднегодовыми числами Вольфа, с учётом полярности цикла; В; = 0,331; Вз = 0,336) 


® для крутящего момента движения Солнца - 2, 3 и 4 года (по максимальным [по модулю] 
корреляциям со среднегодовыми числами Вольфа, с учётом полярности цикла; 
В) = -0,309; Вз = -0,343; В = -0,332) 


В окончательном расчёте использованы данные за 269 лет, соответственно, ошибка коэффици- 
ента корреляции не превышает 0,055; полученные корреляции (и, следовательно, вычислен- 
ные фазовые сдвиги) являются статистически значимыми (Р<0,05). 


2.4. ВЛИЯНИЕ ДВИЖЕНИЯ СОЛНЦА НА ПРОЯВЛЕНИЯ СОЛНЕЧНОЙ АКТИВНОСТИ 


Итак, для расчёта коэффициентов корреляции, способных продемонстрировать влияние грави- 
тационных возмущений (в частности, возникающих при движении Солнца относительно бари- 
центра Солнечной системы) на показатели солнечной активности использованы следующие 
данные (период с 1719 по 2019 год, всего 301 год): 


е значения среднегодовых чисел Вольфа, без учёта полярности цикла; 

е значения среднегодовых чисел Вольфа, с учётом полярности цикла; 

е вычисленные значения «Фактора Е»; 

е вычисленные значения «Фактора Х»; 

е значения отклонений центра Солнца от барицентра Солнечной системы, с лагом 7 лет; 
е значения момента импульса движения Солнца, с лагом / и 8 лет; 

е значения крутящего момента движения Солнца, с лагом 2, Зи 4 года. 


Полученные значения?1 коэффициентов корреляции (с учетом ошибки} приведены в Табл.2. 


Таблица 2 
Корреляции показателей солнечной активности с параметрами движения Солнца 


отклонение 
центра 
Солнца от момент импульса крутящий момент 
барицентра движения Солнца движения Солнца 
Солнечной 
системы 


лаг 7 лет лаг 7лет | лаг 8 лет | лаг 2 года | лаг 3 года | лаг 4 года 


среднегодовые числа Вольфа, -0,020 0,008 0,012 -0,001 0,005 0,007 
без учёта полярности цикла (0,058) (0,058) (0,058) (0,058) (0,058) (0,058) 


среднегодовые числа Вольфа, 0,357 0,370 0,366 | -0,344 | -0,365 | -0,341 
с учётом полярности цикла (0,050) (0,050) (0,050) (0,051) (0,050) (0,051) 
вычисленные значения «Фактора Е» 0,139 0,100 0,128 -0,079 -0,118 -0,137 
(0,057) (0,057) (0,057) (0,057) (0,057) (0,057) 
вычисленные значения «Фактора Х» -0,142 -0,139 -0,165 0,104 0,159 0,186 
(0,056) (0,057) (0,056) (0,057) (0,056) (0,056) 


В скобках приведена величина ошибки соответствующего коэффициента корреляции (М = 301). 





Выделены максимальные [по модулю] коэффициенты корреляции, отражающие несомненное 
влияние параметров движения Солнца на показатели солнечной активности. 
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3. ОБСУЖДЕНИЕ 


Солнечная активность проявляется с определённой периодичностью. Принято считать, что ос- 
новной период составляет около 11 лет-? (впрочем, при ближайшем рассмотрении оказывает- 
ся, что число «11» является «мантрой», а не результатом скрупулезных вычислений). В послед- 
нее время преобладает мнение о 22-летней периодичности солнечной активности, основанное 
на знаниях о смене полярности магнитного поля Солнца и, соответственно, необходимости 
объединять 11-летние «пики» попарно. 


Рассмотрим длительность периодов солнечной активности более детально (Табл.З3). 


Таблица 3 


Продолжительность периодов солнечной активности 
по значениям среднегодовых чисел Вольфа (1700-2019) 


Цикл активности 
(годы) 
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1986 


Отчетливо видно, что реальными наблюдениями не подтверждается преобладание 11-летних 
или 22-летних циклов. Более точно периоды циклов солнечной активности можно вычислить 
по среднемесячным числам Вольфа (Табл.4). 


Таблица 4 


Продолжительность периодов солнечной активности 
по значениям среднемесячных чисел Вольфа (1749-2020) 


Без учёта С учётом 
Цикл активности полярности полярности цикла 


оды, меся 
Боли И - И 


1749.01 — 1755.05 
1755.06 — 1766.06 
1766.07 — 1775.01 
1775.02 — 1783.10 
1783.11 - 1798.06 
1798.07 — 1809.08 
1809.09 - 1822.12 
1823.01 - 1833.06 
1833.07 - 1843.02 
1843.03 - 1855.08 
1855.09 - 1867.01 
1867.02 — 1878.04 
1878.05 - 1889.11 
1889.12 — 1901.11 
1901.12 - 1912.02 
1912.03 - 1923.08 
1923.09 - 1933.08 
1933.09 - 1944.02 
1944.03 - 1954.01 
1954.02 - 1964.07 
1964.08 -— 1975.04 
1975.05 - 1986.06 
1986.07 - 1996.10 
1996.11 - 2008.08 
2008.09 — 2020.02 


< 


17,3 
17,2 
25,7 
23,7 
22,0 
22,5 
23,4 
21,6 
20,3 
20,3 
21,7 
22,1 





ва из 
| 1755.06-1766.06 |110 — 
Корвяьы вт | 8 
| 1775.02-178340 | 87 
ео 
[1798.07 -1809.08 | мл — 
| вл 
[1823.01 -1833.6 | 104 — 
оатьы м | №8 
[1843.03 1855.8 | 124 — 
овтотаты | № 
[1867.02 -1878.4 | м2 — 
Е 
| 1889.42-190414 | 119 — 
ИЕ ИН 
| 1912.03-1923.08 | ид _ 
ато ма | 
[1933.09 - 1944.02 | 104 — 
Кота ма! № 
| 1954.02-1964.07 | 104 _ 
ЕВ ЕО 
| 1975.05-198606 | мл — 
ЕВЕ АЙ 
[ 1996.41-2008.08 | 118 — 
|_2008.09-202002 | 1м | - 


ПЗ 


Очевидно, что и более точные расчёты не подтверждают преобладание циклов продолжитель- 
ностью 11 лет либо 22 года. 


По всей видимости, надо говорить не о близости среднего периода циклов активности Солнца к 
какому-то определённому значению, а о множественном влиянии различных физических фак- 
торов, по-разному проявляющемся в разные эпохи наблюдений. Кроме того, следует соотно- 
сить энергетическую мощность собственных процессов Солнца с мощностью внешних (модули- 
рующих) воздействий. Разумеется, нигде более в Солнечной системе нет электромагнитных по- 
лей, хотя бы сравнимых по мошности с полями самого Солнца. Однако влияние гравитацион- 
ных возмущений может быть заметным, тем более, что оно продолжается в течение несколь- 
ких миллиардов лет. Солнечная система стабильна, орбиты планет за пару-тройку миллиардов 
лет радикально не изменились. Значит, эффект достигнутого резонанса между гравитационны- 
ми возмущениями, приводящими к смещениям Солнца относительно барицентра Солнечной 


системы, а также изменениями момента импульса и крутящего момента может и должен ска- 
зываться на проявлениях солнечной активности. 


Наиболее значимыми среди периодических гравитационных возмущений являются, как пока- 
зано выше, колебания с периодом около 19,8 года (от 2/3 до 3/4 совокупных воздействий). Не 
следует забывать и об эффекте «биений» при наложении нескольких периодических колеба- 
ний. Кроме того, существующая методика обработки наблюдений солнечных пятен и вычисле- 
ния тех или иных значений чисел Вольфа приводит к существенной разнице в динамике соот- 
ветствующих показателей (см. Рис.9). Да, общий характер цикличности наблюдаемых явлений 
очевиден - но при сравнении временных рядов среднегодовых и среднемесячных чисел Воль- 
фа становится заметно не только несоответствие дат пиковых и сглаженных значений, но и не 
всегда чёткое разделение двух циклов, следующих непосредственно друг за другом. 
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Рис.9. Динамика солнечной активности (числа Вольфа): 
черный - среднегодовые значения (1700-2019), красный и синий - среднемесячные (1749-2019); 
чередование цветов отражает смену полярности цикла 


Характерно, что выявленные периоды модулирующих воздействий не совпадают с периодами 
обращения отдельных планет - не какая-то одна планета, или только планеты-гиганты оказы- 
вают влияние на Солнце, но все планеты Солнечной системы, в ТОЙ ИЛИ ИНОЙ степени, в резуль- 
тате сложных гравитационных взаимодействий. 


И неважно, какой физический механизм приводит к спорадическому высвобождению высоко- 
энергетических фрагментов солнечного ядра, постепенно всплывающих к поверхности светила 
и отдающих в виде всплесков электромагнитного излучения свой запас энергии в окружающее 
пространство — важно, что такой механизм существует. И именно на него оказывают модули- 
рующее воздействие гравитационные возмущения, создаваемые всей совокупностью планет 
Солнечной системы. В заключение хочется предположить, что многолетние наблюдения за 
звёздами, при обнаружении периодических (модулированных) выбросов электромагнитной 
энергии, помогут выявить светила, обладающие планетными системами, подобными Солнеч- 
ной. То есть системами, состоящими не только из планет-гигантов, но также имеющие в своём 
составе планеты земной группы на ближних орбитах. 
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ПРИЛОЖЕНИЕ 


Приводится только содержательная часть программ на языке ЯКОБИ; директивы управления 
протоколом, промежуточной печати и проецирования объектов/признаков на фазовые плоско- 
сти не показаны. 


П1. Программа с$ип-уеаг.]ас — вычисление автокорреляций для временных рядов годичных по- 
казателей движения Солнца 


* выполнить * уеаг-аг-согг 902 
* выполнить * уеаг-аг-согг 903 
* выполнить * уеаг-аг-согг 904 
* конец * 


П1.1. Подпрограмма уеаг-аИг-согг.дас — подготовка матрицы данных и вычисление корреляций 
стандартным способом для заданного признака (по годичным данным) 


* очистить память данных * 

* читать данные * с5ип_+{55В_6000у_хл1у 

* ввести объекты с номерами * || [900-908] 

* создать признак * [0] 

* копировать значения признака * ]|[ [$1] [0] 

* исключить фрагмент матрицы данных * [] |0] 
* создать признаки * [1-299] 

* копировать значения признака * || [0] |0] 

* сдвинуть столбцы * ]|[ ]0][ +1 

* исключить объекты с отсутствием данных * || ][ 
* центрировать столбцы * || ][ 

* нормировать столбцы * || || 

* вычислить главные компоненты * ]||[ ][1 

* печатать ковариации(корреляции) * [0] [0-19] 
* * [0] [20-39] 

* * [0] [40-59] 

* * [0] [60-79] 

* * [0] [80-99] 

* * [0] [100-119] 

* * [0] [120-139] 

* * [0] [140-159] 

* * [0] [160-179] 

* * [0] [180-199] 

* * [0] [200-219] 

* * [0] [220-239] 

* * [0] [240-259] 

* * [0] [260-279] 

* * [0] [280-299] 

* конец * 


П2. Программа с$ип-топЕИ.дас — вычисление автокорреляций для временных рядов месячных 
показателей движения Солнца 


* выполнить * топ{-аНг-согг 902 
* выполнить * топ{Й-аКг-согг 903 
* выполнить * топ{-аКг-согг 904 
* конец * 


П2.1. Подпрограмма топЙ-аЦг-согг.дас — подготовка матрицы данных и вычисление корреля- 
ций стандартным способом для заданного признака (по месячным данным) 


* очистить память данных * 

* читать данные * с5ип_+{55В_3000у_х1т 

* ввести объекты с номерами * || ][ 

* создать признаки * [0-599] 

* копировать значения признака * || [$1] [0-599] 
* исключить фрагмент матрицы данных * [] [900-908] 
* сдвинуть столбцы * ]|[ ]0][ +1 

* исключить объекты с отсутствием данных * || ][ 
* центрировать столбцы * || ][ 

* нормировать столбцы * || || 

* вычислить главные компоненты * ]||[ ][1 

* печатать ковариации(корреляции) * ]|[ [0] ][ 

* конец * 


ПЗ. Программа с$ип-М/оН-рпа$е.дас — поиск оптимального значения запаздывания (в годах) ме- 
жду показателями солнечной активности и показателями движения Солнца 


* выполнить * сбип-\М/оЁ-рПазе-согг 902 
* выполнить * с5$ип-\/о|{-рПа$е-согг 903 
* выполнить * сбип-\М/о{-рПа$е-согг 904 
* конец * 


П3.1. Подпрограмма с$ип-М/оМ-рПа$е-согг.дас — подготовка матрицы данных и вычисление кор- 
реляций между показателями стандартным способом 


* очистить память данных * 

* создать признаки * [0,1,2] 

* читать данные * \/оЁ{ 

* ввести объекты с номерами * || [0] 

* копировать значения признака * ][ [0] [1] 

* читать данные * \\/о|{-$12п 

* ввести объекты с номерами * || [0] 

* копировать значения признака * ][ [0] [2] 

* создать признаки * [700,800] 

* читать данные * \\/о{-ЕХ 

* ввести объекты с номерами * ]| [700,800] 

* создать признаки * [900-908] 

* читать данные * с5ип_+55В_1700-2019_х1у 

* ввести объекты с номерами * || [900-908] 

* создать признаки * [100-149] 

* копировать значения признака * ]|[ [$1] [100-149] 
* исключить фрагмент матрицы данных * [] [0,900-908] 
* сдвинуть столбцы * || [101-149] +1 

* исключить объекты с отсутствием данных * || ][ 
* центрировать столбцы * || ][ 

* нормировать столбцы * || || 

* вычислить главные компоненты * ]||[ ][1 

* печатать ковариации(корреляции) * [100-149] [1,2,700,800] ][ 
* конец * 


П4. Программа с$ип-МоЙ-сотр!ех.фас — расчет коэффициентов корреляции между показателя- 
ми солнечной активности и показателями движения Солнца 


* создать признаки * [0,1,2] 

* читать данные * \М/оЁ{ 

* ввести объекты с номерами * || [0] 

* копировать значения признака * ][ [0] [1] 

* читать данные * \\/о|{-$12п 

* ввести объекты с номерами * || [0] 

* копировать значения признака * ][ [0] [2] 

* создать признаки * [700,800] 

* читать данные * \\/о|{-ЕХ 

* ввести объекты с номерами * ]| [700,800] 

* создать признаки * [900-908] 

* читать данные * с5ип_+55В_1700-2019_х1у 

* ввести объекты с номерами * || [900-908] 

* создать признаки * [207] 

* копировать значения признака * ]| [902] [207] 

* сдвинуть столбцы * || [207] +7 

* создать признаки * [307,308] 

* копировать значения признака * ][ [903] [307,308] 

* сдвинуть столбцы * || [307] +7 

* сдвинуть столбцы * || [308] +8 

* создать признаки * [402-404] 

* копировать значения признака * ][ [904] [402-404] 

* сдвинуть столбцы * || [402] +2 

* сдвинуть столбцы * || [403] +3 

* сдвинуть столбцы * ]|[ [404] +4 

* исключить фрагмент матрицы данных * [] [0,900-908] 
* исключить объекты с отсутствием данных * || ][ 

* печатать статистические параметры * ]]| || 

* центрировать столбцы * || ][ 

* нормировать столбцы * || || 

* вычислить главные компоненты * ]||[ ][1 

* печатать дисперсии компонент * || 0.99 

* печатать собственные векторы * ][][ 0.95 

* печатать ковариации(корреляции) * [1,2,700,800] [207,307-308,402-404] || 
* печатать ковариации(корреляции) * [1,2,700,800] [207,307-308,402-404] [1-7] 
* печатать ковариации(корреляции) * [1,2,700,800] [207,307-308,402-404] ]1-7[ 
* конец * 


